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Statistical properties of the antisymmetrized molecular dynamics (AMD) are 
classical in the case of nucleon-emission processes, while they are quantum me- 



\& . chanical for the processes without nucleon emission. In order to understand 



this situation, we first clarify that there coexist mutually opposite two statis- 
tics in the AMD framework: One is the classical statistics of the motion of 



wave packet centroids and the other is the quantum statistics of the motion 
of wave packets which is described by the AMD wave function. We prove 
the classical statistics of wave packet centroids by using the framework of the 
microcanonical ensemble of the nuclear system with realistic effective two- 
nucleon interaction. We show that the relation between the classical statistics 
of wave packet centroids and the quantum statistics of wave packets can be 
obtained by taking into account the effects of the wave packet spread. This 
^ . relation clarifies how the quantum statistics of wave packets emerges from the 

classical statistics of wave packet centroids. It is emphasized that the temper- 
ature of the classical statistics of wave packet centroids is different from the 
temperature of the quantum statistics of wave packets. We then explain that 
the statistical properties of AMD for nucleon-emission processes are classical 
because nucleon-emission processes in AMD are described by the motion of 
wave packet centroids. We further show that when we improve the descrip- 
tion of the nucleon-emission process so as to take into account the momentum 
fluctuation due to the wave packet spread, the AMD statistical properties for 
nucleon-emission processes change drastically into quantum statistics. Our 
study of nucleon-emission processes can be conversely regarded as giving an- 
other kind of proof of the fact that the statistics of wave packets is quantum 
mechanical while that of wave packet centroids is classical. 
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I. INTRODUCTION 



In order to have better description of fragmentation processes in nuclear reactions, sev- 
eral molecular dynamics models with nucleon wave packets have been developed in nuclear 
physics besides the molecular dynamics with point nucleons imported from other fields such 
as molecular physics and solid state physics. An important question about the wave packet 
molecular dynamics is whether its statistical properties can be quantum mechanical unlike 
the point particle molecular dynamics or not. To have proper quantum statistics is highly 
desirable for the molecular dynamics description of particle emission and fragmentation 
processes in relatively low temperature stage of nuclear reactions. 

The purpose of this paper is to clarify in detail the statistical properties of the antisym- 
metrized molecular dynamics (AMD) pl|2[ in general situations including both non-nucleon- 
emission and nucleon-emission processes. Many of our arguments here can be equally applied 
to other kinds of wave packet molecular dynamics. About the statistical properties of AMD, 
there have been proposed two opposite opinions, one proposed by Ohnishi and Randrup || 
who say that the AMD statistics is classical and the other proposed Schnack and Feldmeier 
0] who insist that the AMD statistics is quantum mechanical as well as the statistics of the 
fermionic molecular dynamics (FMD) ||. Both two groups of authors analysed the many- 
fermion system under harmonic oscillator mean field when they deduce their conclusions. 
Under this controversial situation about AMD statistics, the present authors have written 
recently a short paper || in which they have insisted that the AMD statistics is quantum 
mechanical for non-nucleon-emission processes while it is classical for nucleon-emission pro- 
cesses. In the paper the present authors have also shown that an improved treatment of 
nucleon-emission processes which takes into account the momentum fluctuation due to the 
wave packet spread converts drastically the AMD statistics for nucleon-emission processes 
from classical statistics into quantum statistics. 

Our present paper aims to give clear arguments which resolve the above-mentioned con- 
troversial situation between Ref. [Q] and Ref. M and then based on this clarification it aims to 
make the meaning of our previous paper of Ref. more understandable and the arguments 
of it more convincing for the readers. 

One of the basic points of this paper is to elucidate a unique feature of AMD that two 
opposite statistics are coexistent: One is the classical statistics which the centroids of nucleon 
wave packets obey and the other is the quantum statistics which nucleon wave packets obey 
in their motion described by the AMD wave function. A key point to understand this 
coexistent feature of different statistics is to recognize that the temperature of the classical 
statistics of wave packet centroids is different from the temperature which characterizes the 
quantum statistics of wave packets. The fact that the motion of wave packet centroids 
obeys classical statistics was proved by Ohnishi and Randrup || by using the canonical 
ensemble technique in simple cases of many-fermion system under harmonic oscillator mean 
field and also in free gas state. In this paper, we prove this fact by using the microcanonical 
ensemble technique in the case of realistic many-nucleon system interacting with effective 
two-nucleon force. The fact that the motion of wave packets described by the AMD wave 
function obeys the quantum statistics was proved by Schnack and Feldmeier [[| also in the 
case of harmonic oscillator many-fermion system. Their proof was made by calculating 
the occupation probability of the harmonic oscillator single-particle level by performing the 
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long-time average in AMD (or FMD) and then by showing that the occupation probability 
calculated with AMD (or FMD) has just the same value as that given by the quantum 
mechanical canonical ensemble average. In this paper we prove this fact in another way so 
that we can understand the relation between the classical statistics of wave packet centroids 
and the quantum statistics of wave packets. Namely we will show that, when we calculate the 
occupation probability of the harmonic oscillator single-particle level on the basis of the fact 
that wave packet centroids obey classical statistics, the resulting value of the occupation 
probability proves to be almost the same as the value given by the quantum mechanical 
canonical ensemble average. Through this proof we will see that the transition from the 
classical statistics of wave packet centroids to the quantum statistics of wave packets is 
surely due to effects of the wave packet spread. 

In the case of the harmonic oscillator many-fermion system which both Ohnishi-Randrup 
and Schnack-Feldmeier groups treated, the AMD statistics should be regarded as being 
quantum mechanical just in accordance with the opinion of Schnack and Feldmeier. It is 
because the AMD statistics that is relevant for calculated observables is not the statistics 
of wave packet centroids but that of wave packets whose motion is described by the AMD 
wave function. We should note that the calculation of observables is made by the use of the 
AMD wave function. We can say that the AMD statistics in the case of the realistic nuclear 
system is also quantum mechanical so long as we are concerned with dynamical processes 
confined spatially inside the nucleus like the harmonic oscillator many-fermion system. 

However, as we discussed in our previous paper || , the situation changes when we treat 
nucleon-emission processes with AMD. The nucleon-emission process can not be studied with 
the harmonic oscillator many-fermion model nor the free gas model. The AMD statistics 
which is contained in nucleon-emission processes is now not the statistics of wave packets 
but the statistics of the wave packet centroids and hence it is classical. The reason is because 
the nucleon emission process in AMD is governed by the motion of the wave packet centroid. 
In AMD, a nucleon cannot go out from a nucleus if the wave packet centroid cannot go out 
from the nucleus even when the high momentum component of the wave packet can go out 
from the nucleus if treated quantum mechanically. In Ref. || we further showed that we can 
recover the quantum mechanical character of the AMD statistics also for nucleon-emission 
processes by introducing a modified treatment of the nucleon-emission process. An essential 
point of the modification is that we allow a wave packet near the nuclear surface to split 
into high and low momentum components and allow the high momentum component to 
escape from the nucleus. Namely we take into account dynamical effects of the momentum 
fluctuation due to the wave packet spread. In this paper we give a detailed discussion of this 
modification on the basis of our clarified recognition of the coexistence of mutually opposite 
two statistics in the AMD framework. Our study of nucleon-emission processes with the 
ordinary AMD and with the modified AMD can be conversely used as another kind of proof 
of the fact that the statistics of wave packets is quantum mechanical while that of wave 
packet centroids is classical. 

The organization of this paper is as follows. In the next section (Sec. II) we explain 
briefly the AMD formalism. Then in Sec. Ill, we prove that the statistical properties of 
wave packet centroids are classical. The proof is performed by using the microcanonical 
ensemble technique in the case of realistic many-nucleon system interacting with effective 
two-nucleon force. In Sec. IV we prove that the statistical properties of wave packets whose 



3 



motion is described by the AMD wave function are quantum mechanical. The proof is made 
in the case of the harmonic oscillator many-fermion system by calculating numerically the 
occupation probability of the harmonic oscillator single-particle level on the basis of the 
fact that wave packet centroids obey the classical statistics. We will see that the resulting 
value of the occupation probability proves to be almost the same as the value given by the 
quantum mechanical canonical ensemble average. Here we further show that in the case 
of the harmonic oscillator system of distinguishable many particles the proof can be made 
totally analytically. In Sec. V we discuss the statistical properties of AMD for nucleon- 
emission processes. We then discuss a modified treatment of the nucleon-emission process 
which recovers the quantum mechanical character of the AMD statistics also for the nucleon- 
emission process. Summarizing discussions are given in Sec. VI. 



II. BRIEF EXPLANATION OF AMD FORMALISM 

We here explain briefly the AMD formalism. For details the reader is referred to Ref . [|TJ . 
AMD describes the nuclear many-body system by a Slater determinant of Gaussian wave 
packets as 

$(Z) =det[0(r i ,Z fc )x« fc tf)], 
0(r, Z) = (^) 3/4 exp{-i/(r - Z/^f + \l% (1) 

where the complex parameters Z = {Z^} stand for the centroids of wave packets and Xa(j) 
for the spin-isospin wave function with a = p |,p |,n f,n |. The width parameter v is 
common to all nucleons and is time- independent. The time evolution of Z is determined by 
the time-dependent variational principle and the stochastic two-nucleon collision process. 
The AMD equation of motion for Z derived from the time-dependent variational principle 
is 

^ dZ kr &H 

= ^ ZfcT log($(Z)|$(Z)), (2) 

where a,r = x,y, z and the Hamiltonian function TC is the expectation value of the Hamilto- 
nian operator H with the subtraction of the spurious kinetic energies of the center-of-mass 
zero-point oscillations of fragments, 

(#(Z)|TO_«V 

(<$>(Z)\<$>(Z)) 2M 0V F) y ' 

Here A is the total mass number, M is the nucleon mass, Np is the fragment number which 
is a function of Z and Z*, and T is 3h 2 u/2M in principle but is treated usually as a free 
parameter for the adjustment of binding energies. 

For the stochastic two-nucleon process, we introduce physical nucleon coordinates W = 
{W,} as 
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k=l 



Qi*= d{z *. Zk) ^g(Hz)mz)}. (4) 

In the case of light systems where we have less than three nucleons in any spin-isospin 
state, the physical nucleon coordinates W are canonical coordinates. However in other cases 
W are not canonical coordinates. An important property of W is that in the many-body 
phase space of W there exists a Pauli-forbidden region into which any nucleon cannot enter. 
For example, we can easily show that the following region of W is contained in the Pauli- 
forbidden region, 

Ew;'W 3 -<4 (5) 

where Na is the minimum number of the total harmonic oscillator quanta of the A-nucleon 
system. The Pauli blocking of the two-nucleon collision is defined such that the two-nucleon 
collision should avoid entering into the Pauli-forbidden region. 

The AMD equation of motion (Eq. (^)) is not of canonical form, indicating that the co- 
ordinates Z are not canonical coordinates. Although the existence of canonical coordinates 
is known mathematically, we do not know at present the explicit forms of canonical coor- 
dinates in AMD except in the case of light systems mentioned above (just below Eq. (^)). 
However, even if we do not know the explicit forms of canonical coordinates, we can derive 
the explicit form of the Jacobian of the transformation between canonical coordinates and 
the original coordinates Z. If we express a set of canonical coordinates by Y = {Y^}, the 
Jacobian J can be shown to be given by the determinant of the hermitian matrix {Cja^r} 
appearing in the AMD equation of motion (Eq. ([|)), 

a(y,n d(R,p) 

where R = {Rj} and P = {Pj} are real and imaginary parts of Y, respectively, while D = 
{Dj} and K = {Kj} are those of Z: 

In spite of our general ignorance about explicit forms of canonical coordinates, there exists 
one canonical coordinate whose explicit form we know. It is the center of mass coordinate 
Z G , 

1 A 



Dg = ^£% K g = £K,. (8) 
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For the sake of discussions in later sections, we here emphasize the fact that the AMD 
wave function is just the exact solution of the time-dependent Schrodinger equation in the 
case of the harmonic oscillator mean field. Namely, when the Hamiltonian operator H is 
given as 

A 

H = huJ2a] ■ SLj + E , (9) 

j'=i 

where a] and a., (j = 1 ~ A) are harmonic oscillator creation and annihilation operators, 
the exact solution of the Schrodinger equation ih(d/dt)^/ = is given as 

* = e -^det[0(r„e-^zl o) ) Xafc (j)]. (10) 

Here the width parameter v of the Gaussian wave packet </> is v = Mu/2h, and the parameters 
{Z^} stand for arbitrary initial values for Z = {Z^}. Needless to say, {Z& = e~ tuJt zf \ (k = 
1 ~ A)} are the solution of the AMD equation of motion (Eq. (|2|)) with T = 3hu/2M. 
We should note the fact that the AMD wave function which is the exact solution of the 
Schrodinger equation is constructed from the motion of wave packet centroids which make 
just the classical oscillatory motion {Z& = e~ %wt 7if\{k — 1 ~ A)}. This fact constitutes 
the basic background of the character of the AMD statistical properties which we clarify in 
later sections. There we will see that the quantum statistics of AMD is constructed from 
the statistics of wave packet centroids which is classical. 



III. CLASSICAL STATISTICS OF WAVE PACKET CENTROIDS 

We here prove the fact that the motion of wave packet centroids in AMD obeys classical 
statistics. The statistical properties of wave packet centroids are determined by the AMD 
equation of motion. We consider that stochastic two-nucleon collisions are not relevant to 
the character of the statistics but play only the role to hasten the equilibration of the system. 
The AMD equation of motion of Eq. (0) is not of canonical form but it can be converted 
into the canonical form if we use canonical coordinates instead of Z. Hence the dynamics 
of wave packet centroids is classical and the statistical properties of wave packet centroids 
are expected to be classical. However, what inhibits us from concluding straightforwardly in 
this way is the effect of the antisymmetrization contained in the AMD equation of motion. 
As we explained in Sec. II, there exists a Pauli-forbidden region in the phase space of the 
physical coordinates W. Although the physical coordinates W are not always equivalent 
to canonical coordinates, we can guess the existence of a Pauli-forbidden region also in the 
phase space of canonical coordinates. If we have a Pauli-forbidden region, we are not sure 
whether we can simply say or not that the statistics of wave packet centroids is classical. 
Under this situation, we are forced to check explicitly the character of the statistics of wave 
packet centroids. 

In Ref. 0, Ohnishi and Randrup investigated the character of the statistics of wave 
packet centroids by using the canonical ensemble technique in simple cases of many-fermion 
system under one-dimensional harmonic oscillator mean field and also in free gas state. The 
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result they obtained is that the statistics of wave packet centroids is classical. Here in this 
paper, we investigate the the character of the statistics of wave packet centroids by using the 
microcanonical ensemble technique in the case of realistic many-nucleon system interacting 
with effective two-nucleon force. 

In the microcanonical ensemble method, the relation between the energy E and the 
temperature T c is given by 

1 _ d log p(E) 
¥ c ~ dE ' 

P( E )= / (II W - E)5(T> G )5(K G ) (11) 



i (2vrft) J 

I (II (2^3 ) MC]5(H - E)5(V G )6(K G ) 



r- 

A 



3=1 



Here det[C] = det[Cj CTj fe T ] and we have used the relation of Eq. (|) for the Jacobian between 
Z and the canonical coordinates Y . The existence of 5(D G )5(K G ) is because of the fact that 
the center of mass coordinate Z G (Eq. ([8])) is a conserved quantity and makes no contribution 
to the phase volume p(E) of the internal motion of the nucleus. For the sake of the numerical 
evaluation of this relation, we rewrite this relation as explained below We first notice 
the following identity relation 

f) r 6 ^ rfH 
dZdet[C}5(H - E)5(B G )5(K G ) ]T a a 



dEJ / v v ~ dZi 

6A 



where 



A d 3 D d 3 K 
dZ=f[ 3 1 

3=1 



{z a , a = l~QA} = {D„ Kj,j = 1 ~ A}, (13) 

and {a a } are arbitrary functions of Z. The identity relation of Eq. fll2| ) can be easily derived 
by using (d/dE5(H — E))(dH/dz a ) = —d/dz a 5(H — E) and integration by parts. We choose 
{a a } so that they satisfy 

For those {a a } which satisfy Eq. fll4]), we have 

01og<y(D G )<y(K G ) 
2^a Q tt; = 0. (15) 

a os 

By dividing both sides of the identity relation Eq. fll2|) by p(-E) and by using Eq. (|TTD , we 

get 



7 



T = ^ ^ fig) 

where (Q(Z)) stands for the microcanonical ensemble average of the quantity Q(Z) for the 
energy E and is defined as 



dZdet[C]5(H - E)5(B G )5(K G )Q(Z) 

W(z » = 1 m ' (17) 

Under the thermally equilibrated situation of the system, we can equate the microcanon- 
ical ensemble average (Q(Z)) to the long-time average Q{Z); 



(Q(Z)) = Q(Z) } 



1 r*a 



Q(Z) = Km / dtQ(Z(t)). (18) 

t2^00 f 2 — t 1 J tl 



By using this equivalence of (Q(Z)} with Q(Z), we can numerically evaluate Eq. fllq) by 
calculating the time development of Z with the AMD equation of motion for a long time. 
We give in Fig. 1 the calculated relation between the temperature T c and the energy E in the 
system of 12 C. Since the contribution of the term (d/dE)(Y^ a a a d7i/ dz a ) in the denominator 
of Eq. fll6|) is very small, we neglected this term in the calculation. The arbitrary functions 
{a a } were here selected to be 



1 A 



A k=1 
1 A 



a a = z a - — K ka for z a = K ja . (19) 



k=l 



We can easily check that these {a a } surely satisfy Eq. (0). We adopted the Volkov No.l force 
for the effective two-nucleon force. The oscillator width v was chosen to be v = 0.16 fm~ 2 . 
It is easy to see that the relation between the excitation energy E* = E — ^(ground state) 
and the temperature T c is very close to the caloric curve E* /ll = 3T C . Since this means 
that the specific heat remains constant even around zero temperature, the statistics is not 
quantum mechanical. It is well known that the caloric curve E*/A = 3T C is characteristic 
to the classical statistics of the harmonic oscillator many-body system. We can say that 
the reason why the caloric curve in Fig. 1 is closer to E* /ll = 3T C than E*/12 = 3T C is 
because in our present approach by the use of the microcanonical ensemble, the center of 
mass coordinate is not involved in the thermal motion implying that the number of degrees 
of freedom is 3(A — 1) rather than 3A. We thus conclude that the statistics of AMD wave 
packet centroids for the realistic nuclear system is classical. 
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FIG. 1. The relation between the temper- 
ature T c and the energy E for the system of 
12 C. Squares are the calculated results with 
the microcanonical ensemble method. Differ- 
ent values of T c for the same E is due to the 
statistical error. Dotted line shows the relation 
E*/(A - 1) = 3T C with A = 12. 

IV. QUANTUM STATISTICS OF WAVE PACKETS 

The fact that the motion of wave packets described by the AMD wave function obeys 
the quantum statistics was first given an explicit proof by Schnack and Feldmeier || . They 
treated the same model system as Ohnishi and Randrup ||, namely the one-dimensional 
harmonic oscillator many-fermion system. They also treated the harmonic oscillator sys- 
tem of distinguishable particles and found also the quantum mechanical character for the 
statistics of wave packets. Their proof was made by calculating the occupation probability 
of the harmonic oscillator single-particle level by performing the long-time average in AMD 
(or FMD) and then by showing that the occupation probability calculated with AMD (or 
FMD) has just the same value as that given by the quantum mechanical canonical ensemble 
average. Here in this paper, we prove the quantum statistics of wave packets in another 
way so that we can understand its relation with the classical statistics of wave packet cen- 
troids. Our proof is made by treating the same one- dimensional harmonic oscillator systems 
of fermions and also of distinguishable particles. 

Like Schnack and Feldmeier, we investigate the occupation probability of the harmonic 
oscillator single-particle level. In AMD, this quantity which we express as V n (Z) for the 
AMD wave function <&(Z) is given by 

_ mz)\p n mz)) 

p n = c\c n = £ \M^))(M^)l (20) 

where 4> n (x) is the one-dimensional harmonic oscillator eigenfuncion with oscillator quanta 
n. Here the AMD wave function &(Z) should be considered to be a Slater determinant of 
one-dimensional Gaussian wave packet <p(x, Zj) = \j2u/n exp{— v{x — Zjj \fv) 2 + (1/2)Z|} 
without spin-isospin wave function (or belonging to a single spin-isospin state). The investi- 
gation of V n (Z) is made in thermally equilibrated situation of the system. In this section we 
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describe the thermal equilibrium by the canonical ensemble method. It means that we inves- 
tigate the statistical ensemble of AMD wave functions whose centroids Z follow the classical 
canonical statistics. Let T c denote the temperature of the motion of wave packet centroids. 
The average value of V n (Z) under thermal equilibrium, which we express as (V n (Z)) Tc , is 
calculated as 



dZ det[C]V n {Z)e 



-H(Z)/Tc 



dZdet[C}e 



-H(Z)/T C 



{) = WzWz)) ' ( } 

Here the notations, dZ, det[C], and others are obvious one-dimensional analogues of the 
corresponding three-dimensional quantities; for example, 



dZ J ' 

3=1 



2nh 



det[C] = det[^^-log($(Z)|$(Z))]. (22) 

Similarly, the average value of the energy under thermal equilibrium, which we express as 
(H(Z)) Tc , is calculated as 

/ dZ det[C]H(Z)e- n ^ 
(H(Z)) Tc = J— r . (23) 



We have numerically calculated (V n (Z)) Tc and (TL(Z))t c by using the Metropolis sam- 
pling method as in Ref. [^] in the case of the system with fermion number 4 (A = 4). The 
quantity (H(Z))t c is just the quantity which Ohnishi and Randrup calculated in Ref. 0. We 
have ascertained their result in the case of A = 4; namely the calculated value of (TC(Z))t c 
is given almost exactly as 

{H{Z)) Tc = AT c + hu i £{n+ 1 -). (24) 

n=0 A 

We suspect that this relation is analytically exact, since our numerical calculation shows 
that this equality holds with very high accuracy. This result, of course, shows that the 
motion of wave packet centroids obeys the classical statistics. 

Now let us consider the quantity {V n (Z))T c - In classical mechanics, we have no such 
quantity as (V n (Z))T c , because there is no concept of the discrete energy level in classical me- 
chanics. Therefore we have no idea what value the quantity {V n (Z))T c should take in the case 
of classical statistics. One may think that we have the relation (V n (Z))T c oc exp(—nhuj/T c ) 
with some small modification by the Pauli principle. However, we will soon see that the n- 
dependence of the calculated values of (P n (Z))T c is very different from that of exp(— nhiv/T c ). 
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Since the quantity (V n (Z))T c is a quantum mechanical quantity, we should check whether 
the calculated values of (V n (Z))T c are equal or not to what the quantum statistics gives us. 

In the quantum statistics, the average value of the occupation probability of the harmonic 
oscillator single-particle level, which we express as ((Pn))r, is given as 

a (n 1 ...n4P n \n 1 ...n 4 ) a e- E{ni - n4)/T 

n 

T 



(W> " J -'" 



a (n 1 ...r U \n 1 ...r U ) a e- E{ni - n4)/T 



ni...n4 

|m...n 4 ) a = -^det[0 ni (xi)...0 n4 (x 4 )], (25) 

4 1 
E(n 1 ...n 4 ) = Y,( n 3 + n)^- 

Similarly, in the quantum statistics, the average value of the energy, which we express as 
((H)) T, is given as 



a (n\.. .n4 1 ni.. - n i) a E{rii...n 4 )e 



-E(m...n4,)/T 



(W)T = ril "' n4 E V,..n 4 |n 1 ...n 4 re-^-^ " ^ 



n\...ri4 



The necessary and sufficient condition that the system has quantum statistical properties 
is that any of thermally averaged quantities like (V n {Z))T c and (TC(Z))t c are expressed by the 
thermally averaged quantities in quantum canonical ensemble with a certain temperature 
T. What is decisively important here is the recognition that there is no need that the 
temperature T of the quantum statistics is equal to the temperature T c of wave packet 
centroids. The relation between the temperatures T and T c can be obtained for example by 
the requirement that the value of (TC(Z))t c is equal to the value ((H)) t- 

We here summarize as follows what we are going to check now: We first calculate nu- 
merically the relation between T and T c from the relation 

((H)) T = (H(Z)) Tc , (27) 

Then we check whether the relation, 

«P„» T = (V n (Z)) Tc , (28) 

holds or not. In Fig. 2 we show that the numerical calculations assure very well the relation 
of Eq. (]28[) . This result demonstrates clearly that the statistics of the motion of wave packets 
described by the AMD wave function is quantum mechanical. 

Our arguments given above in this section can be made completely analytically in the 
case of distinguishable particles. For distinguishable particles, the AMD wave function is 
simply a product of Gaussian wave packets, 

$(Z) = n0(x„Z,). (29) 
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FIG. 2. Occupation probability of the single-particle level n for the one-dimensional harmonic 
oscillator system with 4 fermions. Squares are the results of (V n (Z))T c calculated by assuming 
the classical canonical statistics of the wave packet centroids. Lines show the quantum mechanical 
relation ((P h ))t- T and T c of each panel are chosen so that they satisfy ((H)) t = (H(Z))t c - 
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The coordinates Z = {Zj} are now canonical coordinates and we have det[C] = 1. The 
calculation of (H(Z)) Tc can be executed analytically, giving the result 

H{Z) = YffiZi + h*", 
i=l 

(H(Z)) Tc = A(T c + huj). (30) 

The calculated result of (TC(Z))t c shows, of course, the classical statistics of the motion of 
wave packet centroids. The calculation of ((H)) t for the quantum canonical ensemble can 
be easily made as shown in ordinary text books of statistical mechanics and we have 

<< g >>^ <*p(Mr)-i + H PD 

From the relation of ((H)) t = (H(Z))t c , we get the relation between T and T c as follows 

exp(hu/T) — 1 

The quantities (V„.(Z))t c and ((P n ))r can be also easily calculated analytically, and we 
have 

A (2j*Z ) n 

-p n (z) = Y. y -^-f L eM-z*z J ), 

^^'^^ d.t/V" (33) 

((P n )) T =A{1- exp{-huj/T)}exv{-nhuj/T). 
In getting the results for (H(Z))t c and (V n (Z))T c , we have used the formula 



i 7* 7 .\ n 

dDjdK .^iL e -<*z* Zje -p Z * Zj p 



i+i 



(34) 



When we use the relation between T and T c of Eq. fl32|), we can easily prove the equality 

(V n (Z)) Tc = ((P n )) T cc exp(-nhu/T). (35) 

Thus we have proved the quantum statistics of the motion of wave packets completely 
analytically in the case of distinguishable particles under the harmonic oscillator mean field. 
This proof in the case of distinguishable particles clearly shows that the transition from 
the classical statistics of wave packet centroids to the quantum statistics of wave packets 
is due to effects of the wave packet spread, because there is no complexity due to the 
ant isy mmet r izat ion . 
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V. AMD STATISTICS FOR THE NUCLEON-EMISSION PROCESS AND 
INCORPORATION OF WAVE PACKET FLUCTUATION INTO AMD 



In previous sections we have explained a very unique feature of AMD that in the AMD 
framework two opposite statistics are coexistent: One is the classical statistics which the 
centroids of nucleon wave packets obey and the other is the quantum statistics which nucleon 
wave packets obey in their motion described by the AMD wave function. We have pointed 
out that a key point to understand this coexistent feature of different statistics is to recognize 
that the temperature of the classical statistics of wave packet centroids is different from the 
temperature of the quantum statistics of wave packets. 

In AMD, the time development of the system is described only through the time devel- 
opment of wave packet centroids which is calculated by the AMD equation of motion and 
the stochastic two-nucleon collision. The effects of the wave packet spread are taken into ac- 
count through the use of the AMD wave function for calculating observables. This ordinary 
AMD framework is expected to work well for the treatment of the internal motion of the 
nucleus which is free from the nucleon emission. It is because, in the case of the model sys- 
tem under the harmonic oscillator mean field where we have no particle emission, the AMD 
wave function is just the exact solution of the time-dependent Schrodinger equation. When 
the ordinary AMD framework is accepted to work well for non-nucleon-emission processes, 
the AMD statistics should be regarded as being not the classical statistics of wave packet 
centroids but the quantum statistics of wave packets, because observables are calculated by 
the use of the AMD wave function which contains effects of the wave packet spread. The 
statistical ensemble of AMD wave functions is expected to describe well quantum statistical 
distributions of single-particle level (and/or momentum) occupation probability. 

However when we look at the description of the nucleon-emission process by AMD, we 
find that the description only through the time development of wave packet centroids is not 
always appropriate. In AMD, a nucleon can be emitted only when the wave packet centroid 
can go out of the residual nucleus. In quantum mechanics, however, even if the wave packet 
centroid can not go out of the nucleus, the high momentum component of the wave packet 
is allowed to go out of the nucleus. Namely, although the wave packet spread should take 
part explicitly in the dynamics of the nucleon emission, the ordinary AMD framework does 
not take account of this role of the wave packet spread. It should be noted that the wave 
packet spread in AMD is large because the energy spread (3/i 2 z//2M) due to the wave packet 
spread is about 10 MeV. The neglect of the nucleon emission due to the high momentum 
component of the wave packet spread will result, for example, in the small cross section of 
emitted nucleons or the slowdown of nucleon emission process in the dynamical stage of the 
reaction. Of course for a high energy nucleon whose wave packet centroid can easily escape 
from the nucleus, the ordinary AMD description has no serious problem. 

According to the above discussion, the emitted nucleons do not reflect the correct mo- 
mentum distribution inside the nucleus. They reflect only the momentum distribution of 
wave packet centroids. It means that the statistical properties of the nucleon-emission pro- 
cess can not be quantum mechanical. 

In order to check the statistical properties of the nucleon-emission process, we investi- 
gated in Ref. || the equilibrated coexistent situation of gas and liquid nucleons. For the 
sake of self-contained discussion, we here recapitulate the method of the investigation in 
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Ref . H . We put A tot nucleons into the potential wall with a large radius and give them the 
total energy E tot . The confinement of nucleons inside the potential wall is made by adding 
to the AMD Hamiltonian of Eq. (|3|) the following term 

IE/CID.-DcmI) (36) 

with 

f(x) = (x — a) 2 9(x — a), 

= Re Zj/yfi, Dcm = 4~ E D i> ( 37 ) 

Atot j = i 

a = 12fm, k = 5MeV/fm 2 . 

Here 6(x) is the step function. When we calculate the time evolution of the system for a long 
time (~ 10 5 fm/c), the system is in a phase equilibrium of liquid nucleons and gas nucleons. 
By liquid nucleons we imply bound nucleons forming a nucleus. The mass number A\ iq of 
liquid nucleons changes with time and we select the moments at which Au q takes a given 
value in order to get the statistical properties of the nucleus with the given mass number. 
We calculate the temperature T g and the long-time average value of the energy Eu q of the 
liquid nucleons. The temperature T g , which should be common to both phases, is calculated 
as the long-time average value of r, where (3/2)r is the kinetic energy (plus the potential 
energy from the wall) per nucleon in the gas phase. Since the gas phase is dilute, the effect 
of the Pauli principle is neglected. When there exist some fragments in the gas phase, they 
are excluded in the calculation of the temperature. In the calculation of the temperature 
T g , i.e. the kinetic energy of gas nucleons, we only use the central value of the wave packet 
Im Zj in the momentum space. It means that we treat gas nucleons as point particles. The 
treatment of gas nucleons as point particles is consistent with the subtraction of the zero- 
point energies of nucleons in the AMD Hamiltonian 7i of Eq. (H). It is to be noted that for 
dilute gas in a large box there is no difference between quantum and classical statistics. The 
single-particle energy spacing in the box with the size 2a is (h 2 /2M)(ir/2a) 2 ~ 0.35 MeV 
which is sufficiently small in our present problem. 

In Ref. |J we showed that the relation between Eu q and T g is very close to the linear 
curve E^ JAii q = 3T g with E\ iq = Eu q — -E/ i(? (ground state). It is displayed again in Fig. 3. 
On the other hand, according to the proof in Sec. Ill, the caloric curve between Eu q and the 
temperature T c of wave packet centroids of liquid nucleons should be very close to the linear 
curve E^ /Au q = 3T C . Thus we obtain the relation of T g = T c . This relation verifies the 
correctness of our conjecture that in the ordinary AMD the emitted nucleons reflect not the 
correct momentum distribution of nucleons but only the momentum distribution of wave 
packet centroids inside the nucleus. 

By using the calculated linear relation between T g and Eu q we can discuss conversely 
as follows. We first regard the gas nucleons as being a thermometer which measures the 
temperature T c of wave packet centroids of liquid nucleons. Namely we regard that the 
temperature T g is just the temperature T c of wave packet centroids of liquid nucleons. It 
is justified because the interaction between the thermometer and the liquid nucleons is 
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FIG. 3. The statistical property of the 
excited nucleus calculated in the phase equi- 
librium of liquid phase and gas phase. Re- 
sults are shown for the usual AMD (squares) 
and for AMD-MF (diamonds). Lines of 
(E*/A)u q = 3T 9 and T 9 2 /(12MeV) are drawn 
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described only by the equation of motion of wave packet centroids. Then, the calculated 
linear relation between T g and Eu q gives us a novel proof of the classical statistics of wave 
packet centroids. This proof is independent of the proof given in Sec. Ill by the use of the 
microcanonical ensemble technique. 

Now we discuss the improvement of the description of the nucleon-emission dynamics so 
that the high momentum component of the wave packet can go out of the nucleus under the 
situation where the wave packet centroid can not go out of the nucleus. For the improvement, 
we have to treat the splitting of the wave packet into high and low momentum parts, one 
going out of the nucleus and the other remaining inside the nucleus. It requires us an 
introduction of a new stochastic process into AMD. There can be various kinds of new 
description of the nucleon emission. From the arguments given above, the essential point 
which any kind of new description should have is that the nucleon emission through the 
splitting of a wave packet is made according at least to the momentum distribution of the 
wave packet. In Ref. 0] we presented a new description named AMD-MF in which we 
apply the new description only when a nucleon is near the nuclear surface and is going to 
be emitted from the nucleus. For more detailed explanation of the AMD-MF the reader is 
referred to Ref. |J. 

When we introduce a new description of the nucleon emission, we have to check whether 
the new description ensures correctly the quantum statistics or not for the nucleon-emission 
process. We discuss the check of the statistics of the AMD-MF also by investigating the equi- 
librated coexistent situation of gas and liquid nucleons. This investigation was made in Ref. 
H entirely in the same way as discussed above for the check in the case of the usual AMD. 
The calculated relation between T g and Eu q by the use of AMD-MF was found to be approx- 
imated by the caloric curve E\ iq j ' A\ iq = T^/(12MeV) with E* iq = E liq — ^^(ground state) 
which is similar to the empirical Fermi gas caloric curve. It is displayed again in Fig. 3. 
This result strongly supports that the temperature T g can be regarded as being equal to the 
temperature T of wave packets of liquid nucleons. We can say that the new description of the 
nucleon emission in AMD-MF ensures proper quantum statistics for the nucleon-emission 
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process. 

Here also we can discuss conversely by using the calculated result of El iq /An g « 
T?/(12MeV). We first regard the gas nucleons as being a thermometer which measures 
the temperature T of wave packets of liquid nucleons. Namely we regard that the tempera- 
ture T g is just the temperature T of wave packets of liquid nucleons. The justification of this 
identification is due to the proper treatment of the nucleon emission in AMD-MF. In AMD- 
MF the nucleon emission is made properly according to the momentum distribution of wave 
packet. Thus we have a good reason to believe that the energy distribution of the emitted 
nucleons reflects adequately the energy distribution of liquid nucleons which is composed of 
the convolution of the momentum spread of wave packet and the momentum distribution 
of wave packet centroids. Then, the calculated caloric curve E* iq /Au q ~ T 2 /(12MeV) gives 
us a novel proof of the quantum statistics of wave packets of liquid nucleons in AMD. This 
proof is not only independent of the proof given in Sec. IV but also more convincing than 
it, because the system treated here is a realistic nuclear system while the systems treated 
in Sec. IV are simple model systems under one-dimensional harmonic oscillator mean field. 
It is to be stressed here that the present proof elucidates most straightforwardly the fact 
that the wave packet spread causes the transition from the classical statistics of wave packet 
centroids to the quantum statistics of wave packets. 

Before closing this section, we comment on the treatment of emitted nucleons as point 
particles in the case of the usual AMD. As mentioned already, this treatment is due to our 
AMD Hamiltonian 7i of Eq. (|3]) where the zero-point energies of nucleons are subtracted. 
The subtraction of the zero-point energy means that the zero-point energy is converted into 
the translational kinetic energy K 2 /2M ,i.e., that the wave packet in momentum space is 
forced to shrink into a plane wave or a point particle. If we do not subtract the zero-point 
energy we can not have low-energy nucleons whose kinetic energy are smaller 

than 3h 2 v/2M. It means that the temperature T g of gas nucleons can not be lower than 
h 2 is/M. Moreover, it means that a nucleon can not be emitted from the nucleus if the 
available energy of a nucleon is smaller than ?>h 2 v/2M. Even when a nucleon can avoid this 
spurious hindrance of nucleon emission, the emission rate is forced to be spuriously reduced. 
Now let us consider the usual AMD with the subtraction of zero-point energies of emitted 
nucleons. Although the subtraction of the zero-point energy favors the nucleon emission 
compared to the case without subtraction, it still underestimates the nucleon emission rate. 
The reason is given just by what we discussed in this section. Namely, in the usual AMD the 
nucleon emission due to the high momentum tail of the nucleon wave packet is neglected. 
Here, in relation with these discussions, we consider the treatment of the width parameter 
v of the wave packet as dynamical variables as in FMD. In such treatment, the result will 
become somewhat between two cases with and without the zero-point energy subtraction. 
Therefore the satisfactory treatment of the nucleon emission will be impossible as long as we 
require the whole wave packet should go out in the process. Finally we comment on the ad 
hoc prescription to put the wave packet spread to the emitted nucleons in the usual AMD 
with the zero-point energy subtraction. This prescription may be regarded as simulating 
the effect of the emission due to the high momentum tail of the wave packet. It surely 
improves the fitting of the calculated momentum distribution of emitted nucleons to data in 
the high momentum region. However this prescription has nothing to do with the emission 
dynamics and for example the slow emission rate remains untouched. If we calculate the 
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gas temperature T g with this ad hoc prescription, the change of T g is simply an increase 
by h 2 v/M and we have no lower temperature than h 2 u/M, showing inadequacy of this 
temperature. 

VI. SUMMARIZING DISCUSSIONS 

We discussed in detail the statistical properties of AMD in general situations includ- 
ing both processes without and with nucleon emissions. The most important point of our 
discussion was to clarify whether the AMD statistics is quantum mechanical or classical. 
The basic background for this point of discussion was as follows: AMD is a kind of wave 
packet molecular dynamics which has been developed in various forms especially in nuclear 
physics in order to cope with quantum mechanical features of the nuclear system. Since 
the point particle molecular dynamics used widely in other fields like molecular physics and 
solid state physics is of classical nature, one tends to think that the wave packet molecular 
dynamics is also of classical nature. Our basic motivation was to show in a clear way that 
this understanding is wrong. We stressed that the strongest suggestion that the wave packet 
molecular dynamics cannot be simply classical comes from the following fact. Namely, the 
exact solution of the time-dependent Schrodinger equation of the many-body harmonic oscil- 
lator system is given by the wave function of the Gaussian wave packet molecular dynamics. 
This fact is well known for the system of distinguishable particles. This fact is, however, 
also true for the system of fermions as can be easily proved. 

In order to discuss the statistical properties of wave packet molecular dynamics, we 
have to clarify the roles of two ingredients of wave packet molecular dynamics, i.e. wave 
packet centroids and wave packet spread. Needless to say, one cannot at all make light 
of the role of wave packet spread in the discussion of statistical properties of wave packet 
molecular dynamics. For example, the momentum distribution of the system can not be 
described without wave packet spread. In general, in wave packet molecular dynamics 
physical quantities are not described simply by using only wave packet centroids but by 
using the wave function including the wave packet spread. 

One of the basic points of this paper was to elucidate the fact that the presence of 
two ingredients gives rise to a very characteristic feature of the AMD statistics. It is the 
coexistence of mutually opposite two statistics in the AMD framework, one being the classical 
statistics which wave packet centroids obey in their motion described by the AMD equation 
of motion and the other being the quantum statistics which wave packets obey in their 
motion described by AMD wave function. The coexistence of the opposite two statistics is 
also true for the wave packet molecular dynamics of distinguishable particles and is a very 
important and novel character of the wave packet molecular dynamics which constitutes a 
decisive difference from the point particle molecular dynamics. 

The fact that the statistics of wave packet centroids in AMD is classical in spite of the 
antisymmetrization was proved in two independent ways in this paper. Both of the proofs 
were different from the proof of Ref. H by Ohnishi and Randrup and were made in the case 
of realistic many-nucleon systems interacting with effective two-nucleon force. One was a 
rather direct proof by the use of the microcanonical ensemble technique. The other was a 
proof by the use of a thermometer made of gas nucleons which measures the temperature of 
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the motion of wave packet centroids inside the nucleus. The reason why the thermometer 
measures the temperature of the wave packet centroids is because the interaction between 
the thermometer and the nucleons inside the nucleus is governed only by the motion of wave 
packet centroids without any effects of wave packet spread. 

The fact that the AMD statistics of wave packets including effects of wave packet spread 
is quantum mechanical was also given two independent proofs in this paper. Both of them 
are different from the proof of Ref. M by Schnack and Feldmeier. One was a proof in the 
case of one-dimensional harmonic oscillator systems of fermions and also of distinguishable 
particles. We calculated the occupation probability of the harmonic oscillator single-particle 
level on the basis of the fact that wave packet centroids obey classical statistics. We then 
showed that the calculated values of the occupation probability were almost the same as the 
values given by the quantum canonical ensemble average. What was decisively important 
in this proof was the recognition of the fact that the temperature of the classical statistics 
of wave packet centroids is different from the temperature which characterizes the quantum 
statistics of wave packets. A merit of this proof is that it clarifies the relation between 
the classical statistics of wave packet centroids and the quantum statistics of wave packets, 
or in other words, it clarifies how the quantum statistics of wave packets emerges from 
the classical statistics of wave packet centroids. The other proof was by the use of a new 
thermometer made of gas nucleons which measures the temperature of the motion of wave 
packets inside the nucleus. The new thermometer interacts with the nucleus now through 
a new mechanism of the nucleon emission from the nucleus. The nucleon emission from 
the nucleus is made not by the ordinary AMD but by the modified AMD named AMD-MF 
which was proposed in Ref. [BJ. In the AMD-MF, the nucleon emission is governed not by 
the motion of the wave packet centroids but by the momentum distribution due to the wave 
packet spread. A merit of this proof is that it clarifies most straightforwardly the role of the 
wave packet spread for the transition from the classical statistics of wave packet centroids 
to the quantum statistics of wave packets. 

The temperature which is relevant to the calculated observables in AMD is of course not 
the temperature T c of wave packet centroids but the temperature T of wave packets. The 
serious mistake of the work of Ref. |j| by Ohnishi and Randrup is that they regarded the 
temperature T c of wave packet centroids as being the temperature relevant to the calculated 
observables. Based on this mistake they insisted a wrong opinion that one should correct 
the AMD framework even for the harmonic oscillator many-body system for which the AMD 
wave function gives the exact solution of the time-dependent Schrodinger equation. 

The conclusion of this paper about the character of the AMD statistics can be said as 
follows: As far as the effects of the wave packet spread are duely taken into account, the 
AMD statistics is quantum mechanical. For processes without nucleon emission, namely 
for the motion confined inside the nucleus, the ordinary AMD takes due account of the 
wave packet spread and the AMD statistics is quantum mechanical without any need of 
essential modification. However, for the nucleon-emission process, the ordinary AMD does 
not duely take account of the wave packet spread. Namely in the ordinary AMD, the emitted 
nucleons reflect not the momentum distribution of nucleons inside the nucleus but that of 
wave packet centroids inside the nucleus. Hence the statistics of the ordinary AMD for the 
nucleon-emission process is classical. In order to recover the quantum statistics also for the 
nucleon-emission process, the description of the nucleon-emission process of the ordinary 
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AMD should be modified so that the nucleoli emission is made according to the momentum 
distribution of nucleons inside the nucleus which is the convolution of the momentum spread 
of the wave packet and the momentum distribution of wave packet centroids. The AMD-MF 
method proposed in Ref. is one of such modified versions of AMD and it was verified that 
that the statistics of the AMD-MF for the nucleon-emission process is actually quantum 
mechanical. 

What is to be noticed in the above discussion is that the reason of the classical statistics 
of the usual AMD for the nucleon-emission process is originating not from any complex 
mechanism but from a simple mechanism of the single-particle motion of emission. Thus we 
have been able to recover the proper quantum statistics for the nucleon-emission process by 
a modification of the single-particle dynamics. In the modification, a proper incorporation 
of the wave packet spread into the dynamical process should be made. Here, however, this 
modification should not be confused with the treatment of the width parameters of wave 
packets as time-dependent dynamical variables. 

The modification of the AMD by the AMD-MF can be generalized so that we can improve 
the description of the cluster emission process and furthermore that of the nucleon-transfer 
process between the projectile and the target in addition to that of the nucleon-emission 
process. We have already formulated such extension and have found that fragmentation 
processes are largely influenced by the extended modification. We will discuss these in other 
papers. 
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